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ABSTRACT 

We present a detailed comparison between far-UV/optical colour Magnitude Di- 
agrams obtained with high-resolution Hubble Space Telescope data and suitable theo- 
retical models for three Galactic Globular Clusters: M3, M13 and M79. 
These systems represents a "classical" example of clusters in the intermediate metal- 
licity regime that, even sharing similar metal content and age, show remarkably dif- 
ferent Horizontal Branch morphologies. As a consequence, the observed differences in 
the colour distributions of Horizontal Branch stars cannot be interpreted in terms of 
either first (metallicity) or a second parameter such as age. 

We investigate here the possible role of variations of initial Helium abundance (Y). 
Thanks to the use of a proper setup of far-UV filters, we are able to put strong 
constraints on the maximum Y (Y max ) values compatible with the data. We find dif- 
ferences AY max ~ 0.02 — 0.04 between the clusters with M13 showing the largest value 
(Ymax ~ 0.30) and M3 the smallest (Y max ~ 0.27). In general we observe that these 
values are correlated with the colour extensions of their Horizontal Branches and with 
the range of the observed Na-0 anti-correlations. 



Key words: Globular clusters: individual (M3, M13, M79); stars: evolution 
zontal Branch; ultraviolet: stars 
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1 INTRODUCTION 

Horizontal branch (HB) stars are the progeny of low-mass 
Red Giant Branch stars (RGB) burning helium in their cores 
and hydrogen in a shell around it (Hoyle & Schwarzschild 
1955). As first noticed by Iben & Rood (1970), the different 
HB star colour distributions observed in old stellar systems, 
is the reflection of the amount of mass lost during the RGB 
phase. 

The scientific community agrees from nearly fifty years 
about the fact that the principal parameter governing the 
shape of HBs in Galactic Globular Clusters (GGCs) is metal- 
licity. The general rule is that metal-rich systems have red 
HBs, while in the metal-poor ones stars are distributed on 
average at higher effective temperatures (bluer colours). Sev- 
eral exceptions have come out during the last decades; re- 
markable cases the cases of NGC6388 and NGC6441 (Rich 
et al. 1997), which despite their metallicity ([Fe/H] ~ -0.6) 
show some of the bluest HBs known among GGCs (Busso et 
al. 2007; Dalessandro et al. 2008). Moreover several clusters, 
sharing similar metal content, reveal different HB morpholo- 
gies, typical cases being the pairs NGC5927 - NGC6388 at 



high metallicities ([Fe/H] ~ -0.4), M3 - M13 at intermedi- 
ate metallicity regime ([Fe/H] ~ —1.5; Ferraro et al. 1997) 

and M15 - M92 at low metallicities {[Fe/H] 2.3). 

These noticeable exceptions have required the introduction 
of a second (Freeman & Norris 1981) and possibly a third 
parameter in order to explain the HB distributions in all 
GGCs. What we can call now the "i-th parameter problem" 
is still a hot topic, as stressed by several authors, we recall 
the reader to Catelan 2009 for a nice review (see also Dotter 
et al. 2010 and Gratton et al. 2010; hereafter D10 and G10 
respectively) Q 

An accurate knowledge of the physical parameters playing 
a role in shaping the HB is extremely important also for an 
appropriate interpretation of distant unresolved stellar pop- 
ulations. In fact it is well known that the HB morphology can 
have a strong impact on the integrated light of stellar pop- 
ulations, affecting colours and line indices (Lee et al. 2002; 

1 Moreover the increase of photometric capabilities has allowed 
to reveal tiny differences and features like gaps along the HB of 
some GGCs (see Ferraro et al 1998 for example), has complicated 
the interpretation. 
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Schiavon et al. 2004; Percival & Salaris 2011; Dalessandro 
et al. 2012). 

Despite the huge efforts made to address this problem, its 
solution is not obvious and still different scenarios are pro- 
posed. One of the reasons that complicates the identifica- 
tion of the mechanisms - other than metallicity - at work in 
shaping the observed luminosity and effective temperature 
distribution of stars along the HB is that there are many 
possible culprits (mass-loss, age, helium abundance ...; see 
Rood 1973 for example) and some of them are not well con- 
strained from theory. 

Age has been identified as the natural global second parame- 
ter by many authors in the past years (Lee et al. 1987, 1988, 
1990; Lee, Demarque & Zinn 1994; Sarajedini & King 1989). 
According to this interpretation older clusters tend to have 
bluer HBs, while younger ones should have on average red- 
der HB morphologies. This scenario appeared in agreement 
with the picture for the Galaxy formation and its early evo- 
lution (Searle & Zinn 1978; Zinn 1985). By means of high 
resolution HST data for a large sample of GGCs, D10 found 
that the existence of outer halo GCs with anomalously red 
HBs fits well the scenario in which age is the second pa- 
rameter. In fact, the behaviour of the 4-5 relatively younger 
clusters in their sample could be reproduced in term of cor- 
relation between age and HB morphology, while the bulk 
of the analyzed targets is peaked around old ages (see ages 
reported by Salaris & Weiss 2002, G10, D10) and doesn't 
show any obvious correlation. Also results by G10 agree on 
the fact that age is the second main parameter driving the 
HB morphology. 

It is also worth noticing that most of these results are based 
on optical CMDs and HB morphology parameters (like the 
well known HBR from Lee et al. 1994), which tend to min- 
imize the importance of blue tails. On the contrary using 
proper combinations of Ultra- Violet (UV) and optical filters 
has an important impact both in term of HB classification 
and comparison with theoretical models. 

Still, age is not able to explain exhaustively the HB mor- 
phology. Detailed cluster to cluster comparisons have shown 
that there are systems with similar iron content and age, 
but remarkably different HB morphologies. A clear example 
is given by the three clusters M3 - M13 - M80, as shown 
by Ferraro et al. (1997, 1998) and at present there is hardly 
a scenario able to give a satisfactory explanation for their 
different morphologies. 

As suggested by Buonanno et al. (1985) and Fusi Pecci et 
al. (1993), age might be one of many and probably the most 
important HB second-parameter, but not the only one. Buo- 
nanno et al. (1995) argued that it is not possible to repro- 
duce the complex HB zoology with a single "global" pa- 
rameter, but more likely we can explain it by thinking of 
a "global combination" of "non-global" quantities and phe- 
nomena related to the formation, chemical and dynamical 
evolution of each cluster. 

The necessity of at least a third parameter transpires also 
from D10 and G10 analyses, in the form of either the lumi- 
nosity cluster density or stellar density (log(p)) - as already 
suggested by Fusi Pecci et al. (1993) - which might corre- 
late with the hot extension of the HBs, or a variation of the 
initial helium abundance (Y), respectively. 
Indeed D'Antona et al. (2005) and Dalessandro et al. (2011 
- hereafter Paperl) have shown that for NGC2808 the main 
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Figure 1. Optical CMDs of the three GGCs considered in our 
analysis. 



parameter that determines the HB morphology is Y. In par- 
ticular in Paperl we have been able to satisfactory reproduce 
the cluster complex HB morphology by assuming three dif- 
ferent sub-populations with He abundances compatible with 
what inferred from the multimodal Main Sequence (MS; Pi- 
otto et al. 2007) and spectroscopic analyses (Bragaglia et al. 
2010; Pasquini et al. 2011). It is worth noticing however that 
only few clusters show evidences of multiple sub-populations 
with different He abundances as inferred from their MS or 
spectroscopic analyses (NGC2808, uu Centauri - Bedin et al. 
2004, NGC6752 - Milone et al. 2010, NGC6397 - Milone et 
al. 2012; NGC1851 - Gratton et al. 2012) 
In order to further investigate the "i-th parameter prob- 
lem", we selected a triplet of GGCs with similar metal- 
licity ([Fe/H] ~ -1.50) and age, namely NGC5272 (M3), 
NGC6205 (M13) and NGC1904 (M79) as templates for a 
comparative analysis of their HBs. These targets are basi- 
cally the same used by Ferraro et al. (1997, 1998), with the 
only exception of M80, which have been replaced by M79 
since recent spectroscopic analyses (Carretta et al. 2009) 
show that it is slightly more metal-poor than the others. As 
shown by Ferraro et al. (1998), M13 has a very extended HB 
with two clear gaps (named Gl and G2 after their analysis); 
while, the HB of M3 covers a much narrower extension in 
effective temperature and shows only a mild indication for 
a gap at the same T e ff as Gl in M13. M79 represents an 
intermediate case between the two. 

As already done for NGC2808 (Paper I) we will use a combi- 
nation of optical and UV photometric images. Details about 
data reduction are presented in Section 2. In Section 3 
we will focus on the long-standing debate about the age- 
differences between these clusters. Comparisons with theo- 
retical models and cluster-to-cluster differences will be dis- 
cussed in Section 4 and 5. 
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2 OBSERVATIONAL DATA SETS AND DATA 
REDUCTION 

The data-set used in the present work consists of Hubble 
Space Telescope (HST) Wide Field Planetary Camera 2 
(WFPC2) images. All targets have been observed in the 
F555W, F336W and F160BW bands, the only exception is 
M3 for which no F160BW data are available, (see Table 1). 
Images and photometry obtained in other bands are also 
available for some of the clusters, however as shown in 
Rood et al. 2008 and Paper I the combination of F555W, 
F3S6W and F160BW seems to be the best filters setup to 
study both RHBs and stars lying at the very hot end of the 
HB. The photometric catalogues of M3 and M13 have been 
already partially presented by Ferraro et al. (1998), while 
the photometry of M79 has been already used by Lanzoni 
et al. (2007). Moreover Far-UV photometry of M79 and 
M13 has been presented by Rood et al. (2008) in order to 
describe in detail the advantages of using the combination 
of F160BW and F555W magnitudes for comparison with 
HB models. 




2 



In all cases the photometric analysis has been performed 
by using ROMAFOT (Buonanno et al. 1983). The reduction 
strategy has been optimized to be particularly sensitive to 
hot stars, as already described by Ferraro et al. (2003). For 
the present analysis attention has been paid to the cross- 
identification of stars detected in the F160BW images. Stars 
observed only in these bands have been force-fitted in the 
other bands (see Paper I). 

The instrumental magnitudes have been calibrated to the 
VEGAMAG photometric system by using zero-points and 
the gain settings reported in Table 5.1 of the WFPC2 data 
handbook Manual Magnitudes thus obtained have been 
also corrected for Charge Transfer Efficiency (CTE) effect 
by means of the prescriptions by Dolphin (2000) and up- 
dated equations listed in the dedicated page of Dolphin's 
web sitqj. 

The resulting colour Magnitude Diagrams (CMDs) for each 
cluster are shown in both the (rriFieoBW, mFimBW — 
rriF5S5w) and (mF336W, m F 336W — m F ss,sw) planes in FigsQ] 
and H 

Star lists have been reported to the absolute coordinate 
system by mean of proper roto-traslations. In particular, 
for each cluster the instrumental coordinates have been re- 
ported to the GSCII astrometric standard system by using 
a large number of primary and secondary standards in com- 
mon, consisting typically of Wide Field ground-based data- 
sets, in the WFPC2 Fields of View (FOVs) (see Lanzoni et 
al. 2007 for details). At the end of the procedure, the typical 
error of the astrometric solution is ~ 0.2". 

2.1 RR Lyrae identification 

The identification of variable stars lying along the HB is a 
crucial step to perform a suitable comparison between ob- 
servations and theoretical models. This is particularly true 
in case of M3, that hosts tens of RR Lyrae stars. 

2 http: / /www. stsci.edu/documents / dhb / web / 

3 http://purcell.as.arizona.edu/ 



Figure 2. Far-UV CMDs of M13 and M79. 

However, none of the data-sets used here is designed for 
variability search, both because of the short time baseline 
covered and number of images acquired. In our CMDs, vari- 
ables are basically sampled at random phase. Thus we de- 
cided to use publicly available catalogues of variable stars in 
order to cross-identify them in our photometric catalogues, 
as already done in Paper I. We retrieved position, magni- 
tude, and classification lists from the "Catalogue of Variable 
Stars in Galactic Globular Clusters (2011)'13 (Clement et al. 
2001). We refer to that paper for all details and references. 
We focused only on RRLyrae stars, neglecting all the 
other different variable types. We first roto-traslated the 
catalogues of variables to our reference frames by using 
CataXcorr, a software developed at the Bologna Observa- 
tory (Montegriffo, private communication). We then looked 
for stars in common with our photometry. All the known 
RRLyrae stars in our FOVs have been recovered in our cat- 
alogues. We found 59 RRLyrae in M3, 7 in M13 and 5 in 
M79 (see Table 1). Their positions in the optical CMDs are 
shown in Figure [T] 



3 METALLICITY AND RELATIVE AGES 

An accurate determination of metallicities and ages is a 
crucial point for HB morphology studies and for the correct 
interpretation of comparative analyses. 

The three GGCs in our sample have a similar iron content: 
according to the recent reassessment of the GC metallicity 
scale by Carretta et al (2009), metallicities are equal 
to [Fe/H]=-1.50 ± 0.05, -1.58 ± 0.04 and -1.58 ± 0.02 
for M3, M13 and M79, respectively. Consistent results 
were also found by Sneden at al. (2004) who obtained 
[Fe/H]=-1.55 ± 0.02 and -1.57 ± 0.07 for M3 and M13 

4 http://www.astro.utoronto.ca/ cclcmcnt/rcad.html 
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Figure 3. Upper panel. Optical CMD of M79 super-imposed to 
M13 after correcting for differences in distance modulus and red- 
dening. Bottom panels. Zoomed view of the TO region. On the 
left, the observed CMDs are shown with two isochrones from 
the BaSTI database with t = 12 Gyr and t = 13.5 Gyr. The 
right panel shows synthetic CMDs obtained from the adopted 
isochrones and the appropriate photometric errors. 

(for homogeneity we report only values obtained from the 
Lick Hamilton spectra). Only in the analysis by Cohen & 
Melendez (2005) M3 and M13 are slightly more metal-rich, 
i.e. [Fe/H]=-1.39 ± 0.02 and [Fe/H]=-1.50 ± 0.01. In 
addition in this case the observed metallicity difference 
(A[Fe/H]~ 0.1) is larger than the uncertainties. 

On the contrary the determination of their ages is a 
more debated argument. In literature there is a quite gen- 
eral consensus that M3 and M79 are old and coeval GGCs. 
There is instead more uncertainty about the age of M13. 
Ferraro et al. (1997) put constraints on the possible differ- 
ence in age between Ml 3 and M3 by comparing their ridge 
mean lines obtained with the optical CMDs used in this 
work, and theoretical models by Dorman et al. (1996). The 
age difference resulted to be smaller than At = 1 — 1.5 Gyr. 
The same conclusions have been reached also by Johnson 
& Bolte (1998) and Rosenberg et al. (1999). According to 
the analysis by Grundahl et al. (1999) M13 may be younger 
than M3 by 0.7 ± 0.2 Gyr, while for Rey et al. (2001) it is 
older by about 2 Gyr. 

Ages reported by other authors for the clusters of our 
analysis agree quite well within the error bars. Salaris & 
Weiss (2002) found a typical difference smaller than At = 1 
Gyr independently of the metallicity scale used. Ages by 
D10 obtained using the isochrone- fitting method are basi- 
cally identical; in fact t=(12.5±0.5) Gyr and t=(13.0±0.5) 
Gyr for M3 and M13 respectively, while M79 is not in their 



Figure 4. As in Fig. [3] but for M3 and M13. 

sample. The same results were obtained by Marin- Franch et 
al. (2009). For G10 instead, M13 is older than M3 and M79 
of about At = 1.75 Gyr. 

It is worth noticing that even moderate differences of the av- 
erage Y among these three clusters do not affect appreciably 
the results of the methods employed to determine their rel- 
ative ages (see, e.g., Marin- Franch et al. 2010). 
We investigated any possible age difference by comparing 
the (mF555w, mF336w — wf555w) CMDs of the three clus- 
ters of our analysis. 

We used M13 as reference and shifted M79 and M3 ac- 
cording to the difference of E(B — V) values and distance 
moduli (see Sections 3.2, 3.3 and 3.4 for details about the 
values adopted) . As apparent from the upper panels of Fig- 
ure [3] and [4] once relative differences in reddening and dis- 
tance are taken into account, the CMDs match well at the 
TO level and along the RGB, consistently with the fact that 
they have almost the same age and share a similar metallic- 
ity. 

In the bottom panels we focus on their Turn Off (TO) re- 
gion. Also from a close inspection, the CMDs have been 
found to nicely overlap at the TO and SGB region. This al- 
ready rules out any significant age difference between M13, 
M3 and M79. However in order to perform a more accurate 
comparison, we also used two isochrones of proper metallic- 
ity and different age from the BaSTI database (Pietrinferni 
et al. 2006). In particular we used as reference the isochrone 
that best fits the observed CMD (t = 12 Gyr) and an addi- 
tional one that is 1.5 Gyr oldei0. As apparent from the left 
bottom panels of Figures [3] and [3] the older isochrone does 

5 Note that the aim of this analysis is note to obtain an estimate 
of the absolute age, but only to give constraints on relative ages 
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not match either the TO or the SGB. 

Starting from these isochrones, and using the estimated pho- 
tometric errors we build two synthetic populations of simi- 
lar size to the observed ones for each comparison. With the 
typical uncertainties of our HST catalogues at the TO mag- 
nitude level, an age difference of 1.5 Gyr would be clearly 
visible and detectable. We then conclude that M13, M3 and 
M79 differ in age by less than 1 — 1.5 Gyr. This result is 
in agreement with estimates by several authors (Johnson & 
Bolte 1998; Rosenberg et al. 1999; Salaris and Weiss 2002; 
Marin-Franch et al. 2009; D10), while it is incompatible with 
the age differences proposed by G10. It is worth noticing that 
an age difference of about 1.5 — 2 Gyr would be compatible 
with the observed data only assuming extremely different 
He abundances (AY ~ 0.1) between these clusters. 



4 THEORETICAL FRAMEWORK 

We have shown that M3, M13 and M79 have similar metal- 
licities and are coeval within typical uncertainties. Therefore 
in this case, the scenario proposed by G10 and D10, identi- 
fying the second parameter with age, is not applicable. We 
therefore investigate here the possible role of variations of 
Helium abundances in the framework of GCs having experi- 
enced multiple formation bursts (D'Ercole et al. 2008; Con- 
roy & Spergel 2011; Valcarce & Catelan 2011) on a short 
timescale (t ^ 100 Myr) in environments enriched by ma- 
terial polluted by AGB stars (Ventura et al. 2002) or fast 
rotating stars (Decressin et al. 2007). As highlighted before, 
Helium has been proposed also by G10 to be the possible 
HB third parameter. 

As in Paper I, we have compared our photometric data with 
a suitable set of BaST0 a-enhanced HB tracks (Pietrin- 
ferni et al. 2006). We have considered theoretical models 
for [Fe/H] = — 1.62, which is the closest value available in 
the BASTI database to the clusters spectroscopically deter- 
mined [Fe/H]. 

In our analysis we take advantage of all initial He- 
abundances available in the BaSTI database, i.e. Y=0.246, 
0.300, 0.350 and 0.400. Bolometric corrections and extinc- 
tion effects have been calculated as in Paper I. 

Recently, Gratton et al. (2011) have found spectroscop- 
ical evidence of a moderate Na-O anti-correlation in stars 
belonging to the HB of NGC2808, at T cft below ~12000, 
that marks the onset of radiative levitation effects (see dis- 
cussion below). A similar result has been found for M 4 by 
Marino et al. (2011), the general trend being that the bluer 
HB stars are on average more Na-rich and O-poor. As a con- 
sequence of these empirical results, we verified by perform- 
ing a number of tests (see Appendix A) how adequate is the 
use of 'standard' a-enhanced spectra to determine bolomet- 
ric corrections in this temperature regime. We found that 
results obtained by using standard a-enhanced bolometric 
corrections, are unaffected by the presence of the observed 
CNONa anti-correlations. 

Regarding the effect of radiative levitation, spectro- 
scopic investigations by Behr (2003), Moehler et al. (2003), 



6 Available at the following URL: 

http: / /www. oa-teramo.inaf.it/BASTI 



Fabbian et al. (2005) have disclosed that - similar to the case 
of NGC2808 - in these three clusters HB stars hotter than 
~ 12000 K display surface abundances of metal species like 
Fe, Cr and Ti around solar or higher, due to radiative levita- 
tion in the atmosphere. Also, the surface He abundances are 
depleted by a factor of ten or more as a consequence of grav- 
itational settling. Following the approach used in Paper I, 
we mimic this effect using 'standard' a-enhanced HB stellar 
models by applying bolometric corrections appropriate for 
[Fe/H]=0.0 (and scaled-solar mixture) when T e fi is above 
12000 K. This is of course a crude approximation, which we 
are forced to adopt because of the lack of extended grids 
of both HB stellar evolution and atmosphere models with 
a large range of chemical compositions, that include consis- 
tently the effect of radiative levitation. The underlying as- 
sumptions are that (%) the HB model evolutionary lifetimes, 
(it) the evolution of T c h and (Hi) bolometric luminosity, and 
(in) the structure of the corresponding model atmospheres 
are not affected by the radiative levitation, and that the rel- 
evant bolometric corrections are determined mainly by the 
enhanced metals. With these assumptions the effect of ra- 
diative levitation on bolometric corrections makes rriFieoBW 
fainter by about 0.15 and m_F336W and rriFhhhw brighter by 
about 0.1 and 0.05 respectively. Paper I has shown how this 
approximation allowed us to recover the three different He- 
abundances inferred by the multimodal MS, along the blue 
HB of NGC2808. We are therefore confident that a similar 
procedure will be adequate to constrain the He distribution 
along the HB of these three clusters, that we are going to 
discuss separately in the following section. 

4.1 The synthetic HB simulations 

The theoretical analysis performed in this paper is based 
on synthetic HB calculations, performed as described in Pa- 
per I. The analysis of NGC2808 was simplified by the fact 
that the initial He content of the cluster subpopulations is 
somewhat 'quantized', as shown by the observed trimodal 
MS (Piotto et al. 2007). In case of M13, M3 and M79 there 
is no apparent 'quantization' of the MS, and we need to con- 
sider the possibility of a continuous distribution of Y along 
the HB. 

We stress here that in these cases the modeling of the Y 
variations along the HB becomes more uncertain since cal- 
culations require more assumptions. Therefore the main goal 
of this work is not to attempt an exact description of the 
Y distribution, that would be affected by our assumptions, 
but to understand whether helium plays a role in shaping 
these very different HBs. UV filters are much better suited 
than optical ones to address this issue. 

Our synthetic HB calculations require the specification 
of at least 4 parameters, plus the cluster age, that we assume 
to be equal to 12 Gyr (see Section 3). These four parame- 
ters are the minimum value of Y (Y m ; n ), the range of He 
abundances (AY), the mean value of the mass lost along 
the RGB AM - that in first instance (and for simplicity) we 
assume to be the same for each Y - and the spread around 
this mean value (cr (AM)). Throughout our analysis we will 
assume that Y varies according to either a uniform prob- 
ability distribution, from Y m i n to Y m i n + AY, or a Gaus- 
sian one. In this latter case AY is the la spread around 
a prescribed mean value. Also for the RGB total mass loss 
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we will assume a Gaussian distribution around the mean 
value AM or a uniform distribution with minimum value 
equal to AM. The idea behind this type of simulations (see, 
e.g., Caloi & D'Antona 2005) is that the colour extension 
of the HB is driven mainly by the variation of Y rather 
than mass loss efficiency. A higher initial Y implies a lower 
mass at the tip of the RGB (for a fixed cluster age) hence, 
for a fixed value of AM, a smaller mass along the HB and 
a bluer colour (the He-core mass is affected by the varia- 
tion of Y to a smaller extent). Our synthetic HB code first 
draws randomly a value of Y, and determines the initial 
mass of the star at the RGB tip (Mtrgb) from interpo- 
lation among the BaSTI isochrones of the prescribed age. 
The mass of the corresponding object evolving along the 
HB (Mhb) is then given by Mhb = Mtrgb — AM plus a 
Gaussian random perturbation <t(AM) (or a uniform prob- 
ability between AM and AM+a(AM)). The WFPC2 mag- 
nitudes of the synthetic star are then determined according 
to its position along the HB track with appropriate mass 
and Y- obtained by interpolation among the available set of 
HB tracks - after an evolutionary time t has been randomly 
extracted. We determined t assuming that stars reach the 
ZAHB at a constant rate. We employed a flat probability 
distribution ranging from zero to tHB, where tHB denotes 
the time spent from the ZAHB to the He-burning shell ig- 
nition along the early AGB. The value of tHB is set by the 
mass with the longest lifetime (the lowest masses for a given 
Y and Z). This implies that for some synthetic object the 
randomly selected value of t will be longer than its tHB or, in 
other words, that they have already evolved to the next evo- 
lutionary stages. Finally, the derived synthetic magnitudes 
are perturbed with a Gaussian la error determined from the 
data reduction procedure. 

We detail now the effects on these results of possible 
variations of age and chemical composition within the uncer- 
tainties described in Section 3. The mass distribution along 
the HB (as derived in our analysis) does not depend on the 
cluster age, but changing the cluster age changes the mass 
loss necessary to reproduce the HB mass distribution. An 
increase of age of 1 Gyr leads to a decrease of the total mass 
loss along the RGB by ~ O.O2M0 (because of the decreased 
evolving mass). 

One can simulate the effect of a variation of [Fe/H], [a/Fe] 
and/or their sum by considering that at these metallici- 
ties scaled solar and a-enhanced models (also for the HB) 
with the same total global metallicity [M/H] are coinci- 
dent (Salaris et al. 1993). We have therefore considered a 
change A[M/H] = +0.3 dex that corresponds to changes 
of [Fe/H] or [a/Fe] or [Fe/H] + [a/Fe]. We have two effects 
here. The first one is on the derived mass distribution (at 
fixed Y) along the HB. For the redder part of the HB 
we use the (mF336w,mF336w — tiifsssw) colour. At fixed 
colour in the horizontal part of the (mF336W, m,F336w — 
mF555w) CMD, an increase A[M/H]=0.3 decreases the de- 
rived value of Mhb by ~0.04 Mq. For the blue tails we 
use the (m F iu)BW, m F iaoBW - m F z,bz,w) CMD. At fixed 
(rriFiwBW — TnF555\v) an increase A[M/H]=0.3 decreases 
the derived value of Mhb by ~ 0.01 Mq . This means that 
if a larger mass loss is needed to explain the blue HB stars 
(compared to the mass loss needed to explain the redder 
population), an increase of [M/H] tends to reduce this mass 
loss difference. The second effect is on the RGB mass (hence 
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Figure 5. Observed (left panel) vs synthetic (right panel) 
( m F336W , m F336W ~ m F555w) CMD of M13 HB. ZAHB se- 
quences for Y=0.246 (solid line), 0.300 (dashed line) and 0.350 
(dot-dashed line) are also displayed. The rectangular box marks 
the RHB region (see text for details). 



on the absolute value of the RGB mass loss). An increase 
A[M/H]=0.3 causes an increase of M RG b by ~ 0.025 Mq. 
The distribution of Y along the observed HBs is basically 
unaffected by the precise age and [M/H], because it is based 
on the variation of magnitude (in m,F336W and rriFieoBW , 
depending on the star colour) with Y, a differential property 
that is unaffected by the exact model metallicity and age. 



5 COMPARISON BETWEEN THEORY AND 
OBSERVATIONS 

In this section we will discuss separately the theoretical anal- 
ysis of our three selected clusters. We stress again that al- 
lowing Y to vary within individual clusters, increases the 
number of free parameters needed to match the observed HB 
star distribution in the CMD. Leaving aside for a moment 
the open problem of the RGB mass loss, the distribution of 

Y among stars fed onto the HB will depend on the convolu- 
tion of the unknown (at least for these three clusters) initial 

Y distribution among MS stars, with the Y-dependent evo- 
lutionary timescales along MS and RGB predicted by stellar 
models. As discussed before, in absence of constraints from 
the MS, we assumed here as a first order approximation 
either Gaussian or uniform probability distribution for Y 
among HB stars. As a consequence, the derived distribution 
of Y and mass range (at fixed Y) are most probably just a 
'reasonable' approximation to the real ones. A much more 
solid result is the total range of Y covered by the stars in 
each individual cluster, because - thanks to our UV CMDs 
- this is strongly constrained by the observed stellar magni- 
tudes rriFi&QBW (when available). 
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5.1 M13 

Spectroscopic analyses of this cluster reveal a well devel- 
oped Na-O anti-correlation (Sneden et al. 2004, Carretta et 
al. 2009), comparable to the case of NGC2808, but photo- 
metric analyses currently available do not show evidences 
for a multimodal MS. Here we analyze our HB photometry 
with - as stated before - the main aim of constraining the 
range of initial He values of the cluster stars. The use of 
the F160BW filter - not considered in the investigations by 
Caloi & D'Antona (2005) and D'Antona & Caloi (2008) - 
will set strong constraints on the maximum value of Y. 

The first step of our analysis consists in defining the 
highest possible Y value necessary to fit the observed HB. 
To this purpose we calculated a synthetic HB with this pre- 
liminary set of parameters (uniform distributions for Y and 
Gaussian for M ): Y min =0.246, AY=0.104, AM=0.21M o , 
(t(AM)=O.O25M0 . Even though this selection of parame- 
ters is not optimized for a perfect fit to the observed stellar 
distribution along the HB, it will suffice to put a first gen- 
eral constraint to the maximum possible range of initial Y. 
The number of stars in this simulation is higher than the 
observed sample, to reduce the effect of statistical number 
fluctuations on the synthetic HB. 

The parameter AY has been arbitrarily fixed to reach 
the maximum He mass fraction Y=0.35. Once Y m i n is cho- 
sen, the mean RGB total mass loss AM is constrained by 
matching the red boundary of the observed HB (the subpop- 
ulation with the lowest Y has the largest HB mass at fixed 
age and AM). The value of <t(AM) then introduces a general 
dispersion in the relationship between HB colour and local 
He abundance. The smaller the value of this parameter, the 
narrower the range of Y in a given colour bin along the HB. 
We notice that a very small value of cr(AM) (a(AM)=0.005 
for example), would produce a sloped HB (brighter towards 
bluer colours because of a very tight colour-Y relation- 
ship) in the (m,F336w, mF336w — mFS55w) CMD at colour 
(mF33ew—mF555w)> — 0.5, instead of the observed roughly 
horizontal structure (see Figure 0. 

As done in Paper I, we determined the cluster dis- 
tance modulus by fitting the synthetic data to the ob- 
served Red HB (RHB) (defined here as the region of the 
HB contained in the square box plotted in Figure [5} in the 
(mF336W, niF336w — mF555w) plane. Notice how the RHB 
is roughly horizontal in this CMD, and allows us to esti- 
mate straightforwardly the cluster distance by comparing 
the synthetic data to the observed number distribution of 
objects in m,F336w, assuming a reddening E(B— V)=0.02 
(see for example Ferraro et al. 1999). A fit to both the 
lower envelope - defined as the faintest magnitude bin where 
the star counts are at least 2a above zero - of the ob- 
served binned distribution in m,F336w, and the mean value 
of the m,F336W RHB magnitudes (< m,F336w >=15.26) 
provides (m - M) = 14.32±0.o£|, in fair agreement within 
the uncertainties with Ferraro et al. (1999; (m — M)o = 
14.43). The left and right panels of Figure [5] show the ob- 
served and simulated (rriF336W, mF336W — niF555w) CMD, 
respectively. We have also over-imposed three ZAHBs for 
Y=0.246 (solid line), 0.300 (long dashed line) and 0.350 
(dot-dashed line). Figure |S] is equivalent to Figure [S] but for 

7 The error bar is equal to half the width of the adopted bin size 
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lil F160BW 1I2 F555W lil F160RW lil F555W 

Figure 6. As in Figure \E\ but for the {mpieoBW > m-FWOBW — 
m F555w) CMD (see text for details). In the right panel we also 
show two selected HB tracks corresponding to M ~ 0.57Mq and 
M ~ O.5OM for each ZAHB. 

the (m F i60BW,niFiaoBW - rriFasaw) CMD. Whilst in the 
(mF336W, wf555w) CMD ZAHB sequences and synthetic 
stars with differing initial Y are largely degenerate along 
the bluest HB tail as a consequence of the huge increase of 
the bolometric correction in the optical filters, this same tail 
appears as a roughly horizontal sequence at (rriFieoBW — 
mF555\v)< —1-0 in the (mFieoBW, rriFieoBW — rn,F555w) 
CMD. Theoretical ZAHB sequences display here different 
shapes at varying Y (see also Pietrinferni et al. 2006). It is 
clear from a comparison of the observed and synthetic HB 
that Y=0.35 for the bluest stars is too high. The lower en- 
velope of the observed stellar distribution lies around the 
Y=0.30 ZAHB. This implies that the maximum Y (Y max ) 
is -0.30. 

As a conclusion, in the hypothesis that the colour ex- 
tension of the HB is due to a range of initial Y values, 
Y m in=0.246 at the red end of the observed HB implies 
Ymax ~0.30 at the blue end, i.e. the initial He-abundance 
spans a range AY ~0.05. This constraint on Y max is very 
general; the critical discriminating factor is provided by the 
rriFieoBW magnitude of the ZAHB at the blue end of the 
HB. 

One can of course try a more accurate fit to the stellar 
distribution along the observed HB by comparing theoreti- 
cal and observed histograms of star counts as a function of 
colour and magnitude. In this case a fine tuned calibration 
of the set of free parameters is necessary. In the assumption 
(for simplicity) of either uniform or Gaussian probability 
distributions for the mass loss and initial Y, one needs to 
split the simulation into three different sections, because no 
single continuous distribution of Y and mass allows to match 
the observed CMDs. As a general method - applied also to 
M3 and M79 -, for each choice of the free parameters we have 
produced 100 synthetic distributions, each of them with the 
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Figure 7. As in Figure \E\ but for one realization of the best fit 
synthetic HB model discussed in the text. 

same number of objects as in the observational sample. From 
this ensemble of simulations we have determined the mean 
value and associated error for each colour or magnitude bin 
used in their comparison with data. 

The star counts along the RHB on the red side 
of the gap at (m,F336w — rn,F555w)^ — 0.5 are best 
fit with Y min =0.246, AY=0.01 and uniform probabil- 
ity distribution for Y, AM=O.21M0, with a Gaussian 
spread <t(AM)=O.O1M0. The resulting distance modulus is 
(m — M)o = 14.30±0.05, quite similar to the value obtained 
from the simulation in Figure [5] 

The second section of the simulation covers the re- 
gion from the blue boundary of the gap to a colour 
(rriFieow — ?riF555w)~ — 3.0 along the HB blue tail. In 
this case one needs a Gaussian distribution of Y, with mean 
value < Y >=0.285 and la spread A(Y)=0.012, mass loss 
AM=0.235 M with a Gaussian spread cr(AM)=O.OlM . 
The last section of the simulation covers the extreme tail of 
the HB, that is best reproduced with < Y >=0.300 with 
a Gaussian la spread A(Y)=0.003, and AM=0.266 M 
with a Gaussian spread <t(AM)=O.OO2M . Figure[7]and Fig- 
ure [5] display a comparison between the observed HB in the 

(mF336W, W-F336W — TriF555w) and (m,F160BW, tufwobw — 

mF555w) CMDs, and one realization of the simulations with 
the parameter choice described before. From the parame- 
ters of this simulations, it is clear that a good match to the 
observed star counts can be obtained only by relaxing the 
assumption of constant AM for all initial He contents. In 
particular, it appears necessary to slightly increase on aver- 
age AM when Y increases, as already highlighted in Paper I 
for the case of NGC2808. However it is important to note 
that, the maximum Y range is fully constrained indepen- 
dently of the assumptions on mass loss. 

Note that the slight increase of AM with Y could be 
in principle due to an age variation at increasing Y. In this 
scenario the variation of AM corresponds to a variation of 



: m 13 

13.5 - 




-4 -2 0-4-2 

m F160RW~ m F555W m F160BW~ lrl F555W 

Figure 8. As in Figure 151 but for one realization of the best fit 
synthetic HB model discussed in the text. 

Mtrgb with Y, keeping the total RGB mass loss unchanged. 
However, these mass differences imply an increase of age 
with increasing Y (lower mass evolving along the RGB) of 
the order of ~ 4Gi/r that is hard to justify. In fact such 
a large age difference would be detectable from the cluster 
TO and sub-giant branch regions, at odds with the empirical 
evidence (Figure [3] and Figure [4j. 

The mean He abundance along the HB of this simu- 
lation is < Y >=0.270, with a la spread cr(Y)=0.02. No- 
tice that Salaris et al. (2004) determined for this cluster 
Y=0.285±0.024 from the R-parameter, consistent with the 
results of the simulation. 

We finally address the question of what happens if Y m i n 
in RHB stars is increased. As a test, we have first consid- 
ered a constant Y=0.280 for the RHB - as proposed by 
Caloi & D'Antona (2005) and D'Antona & Caloi (2008) 
- and produced the corresponding synthetic population by 
tuning appropriately the value of the mass range along the 
RHB. The distance modulus derived from the rriF336w mag- 
nitudes of RHB stars is 0.14 mag larger than the case of 
our simulation with Y m i n =0.246. Given that the rriFimBW 
magnitude along the extreme blue tail of the observed HB 
increases with increasing Y (see Figure [6| , this higher dis- 
tance modulus would imply Y~0.25 along the extreme HB 
((rriFieoBW — mFttsw) < —3.5), i.e. a value lower than along 
the RHB, a result justifiable only assuming a strong anti- 
correlation between mass loss and Y. 

This test also implies that in principle, for a fine tuned Y m i n 
intermediate between 0.246 and 0.280, the whole HB of M13 
can be also modeled - at least from the point of view of 
the magnitude levels in rriF336w and rriFieoBW - employing 
a single value of Y. We find that when Y m i n ~0.265, the 

(rn,F336W, TUF336W — flF555w) and (mFWOBW ,TUF160BW — 

tuf555w) observed CMDs and star counts can be reproduced 
with a constant He abundance, for a total mass spread of 
~0.14 Mq. Even if we cannot discard this scenario on the 
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Figure 9. Observed (left panel) vs one realization of the best 
fit synthetic (right panel) (m F336w , m F336w - m F555W ) CMD 
of M3. The ZAHB sequences for Y=0.246 and Y=0.300 are also 
displayed. The rectangular box marks the region considered for 
the fit (see text for details). 



basis of the HB analysis, we consider it very unlikely. In 
fact it would mean that M13 lost completely its first gen- 
eration of standard He stars, while preserving most of the 
He-enriched population. This would be hardly compatible 
with results from N-body simulations (see for example De- 
cressin, Baumgardt & Kroupa 2008). 



5.2 M3 

Spectroscopic analyses of M3 RGB stars reveal a Na-0 anti- 
correlation that is reduced in extension compared to M13 
(see, e.g., Sneden et al. 2004, Johnson et al. 2005 and ref- 
erences therein). Also in this case, the currently available 
photometric data do not show evidence of a multimodal 
MS. Also, as well known, the HB of M3 displays a much 
shorter colour extension compared to M13. Here we dis- 
cuss its properties only in the (m,F336W , rriF336w — wifsssw) 
CMD, because no F160BW images are available for this 
cluster. Given the way in which HB stars are distributed in 
all CMDs involving F255W , it doesn't give any additional 
information respect to the optical plane. Also given the rel- 
atively short extension of the HB of M3, we decided not to 
use the shorter wavelength F255W filter. 
We considered in our analysis all objects within the rect- 
angular box in Figure [9] Our synthetic calculations show 
that all these stars have T e ff below the threshold for the on- 
set of radiative levitation. There are three stars outside the 
box, that may be genuine HB objects, and appear spread 
along the bluest tail of the cluster HB, that is vertical also 
in this CMD. Given their extremely small number compared 
to the total sample of HB stars, and the fact that their loca- 
tion does not provide any further constraint on the cluster 



Figure 10. Mass and Y distribution as a function of the 
( m F336W ~ m F555w) colour along the HB of M3 (grey crosses) 
and M13 (black dots). 

He distribution, we do not consider them in the theoretical 
analysis. 

Following the same philosophy as for the case of M13, in 
the assumption of either uniform or Gaussian distributions 
for the mass loss and initial Y, observed star counts as a 
function of the m,F33ew magnitude and (rriF336W — mF555w) 
colour are best reproduced with Y m i n =0.246, AY=0.02 and 
uniform probability distribution, a minimum RGB total 
mass loss AM=0.122Mq with a range (again uniform prob- 
ability distribution) of values <r(AM)=9.0 x (Y-0.246) M , 
where Y is the actual initial He mass fraction of the syn- 
thetic star. It is obvious that this dependence of a(AM) 
on Y - not imposed a priori but constrained by the fit to 
star counts as a function of magnitude and colour - im- 
plies an increase of the mass loss spread with increasing Y, 
that is required by the observed shape of the CMD in the 
{rriF336W, rriF336W — rn,F555w) plane. 

In fact as shown in Figure the reddest part of the 
ZAHB describes a steeply sloped sequence, before turning 
horizontal at rriF336W ~ 16.0. To match the observed star 
counts we need to increase the mass loss spread with increas- 
ing Y. This is necessary to reproduce both the steep increase 
of m,F336w with decreasing colour when (rriF336W — mF555w) 
is larger than ~0.2 - e.g., the lower He subpopulations have 
to cover a relatively small mass range - and the almost con- 
stant m,F33ew when the colour (rriF336w — tifsssw) is bluer 
that this limit - e.g., the higher He subpopulations have 
to cover an increasingly larger mass range. The mean He 
abundance obtained from the synthetic HB is < Y >=0.256, 
with a ler spread cr(Y)=0.006, consistent with the value 
Y=0.253±0.013 obtained from the R-parameter by Cassisi 
et al. (2003). In the synthetic HB calculations we have dis- 
regarded all objects lying in the RRLyrae instability strip, 
whose red and blue edges have been fixed at, respectively, 
logT c ff=3.80 and logT c ff=3.88 (see, e.g., the discussion in 
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Figure 11. Observed (left panel) vs one realization of the best 
fit synthetic (right panel) (m F336w ,m F336w - m F555W ) CMD 
of M79 HB. The box marks the region used for the estimate of 
the cluster distance (see text for details). 



Figure 12. As in Figure [TT1 but for the {mpieoBW > m Fl6QBW — 
m F555\v) CMD (see text for details). 



D'Antona & Caloi 2008), which define a very narrow colour 
width of A(rriF336w — rriF555w) < 0.05 at fixed mF 33 @w ap- 
proximatively in the range 0.35 < (mF3mw — rn,F555w) < 
0.50 . 

The distance modulus obtained as described in the case 
of M13 is (m - M) = 15.00±0.04 with standard reddening 
value E(B— V)=0.01 which is fully in agreement with pre- 
vious estimates (Ferraro et al. 1999; Harris 1996 - cata- 
logue version 2010). In this case the maximum value of Y 
(Y max = 0.266) is essentially constrained by the brightest 
limit of the stellar distribution at (rriF336W — rriF555w) ~0.2. 

Figure [10] displays stellar mass and Y distributions as a 
function of the (m,F336w — mF555w) colour in one represen- 
tative synthetic realization of the cluster HB, compared to 
the corresponding best-fit composite simulation of M13. Not 
only the range of He values is obviously much larger in M13, 
but also the shape of the Y trend with colour is completely 
different in these two clusters. The mass distribution forms 
a very smooth continuous sequence, with the maximum and 
minimum masses along the HB distribution in M3 being ap- 
preciably larger than the corresponding values in M13. 

A detailed comparison with recent results by Catelan 
et al (2009) and Caloi & D'Antona (2008) is reported in 
Appendix B. 

5.3 M79 

Spectroscopic analyses of M79 RGB stars disclose an exten- 
sion of the Na-O anti-correlation that is intermediate be- 
tween that of M13 and M3 (see Sneden et al. 2004, Carretta 
et al. 2010), and photometric analyses do not reveal evi- 
dences of a multimodal MS. Also the colour extension of 
its HB is intermediate between M13 and M3. Following the 
same approach as for M13 and M3, we find that observed 



star counts as a function of magnitudes and colours are best 
reproduced with Y m i n =0.246, AY=0.035 and uniform prob- 
ability distribution, a RGB total mass loss (at a given Y) 
ranging from AM=O.11M to AM=0.27M Q , with a uniform 
probability distribution. One realization of the best fit sim- 
ulation is displayed in Figures [TT] and 1121 

The distance modulus has been determined from the 
fit to the stars contained within the box displayed in 
Figure [TT] (denoted here as RHB) as described for M13. 
Assuming E(B— V)=0.01 (Harris 1996) we have obtained 
(m — M)o=15.64±0.05, values well in agreement with both 
Ferraro et al. (1999) and Harris (1996). For the RHB 
stars we determine an average < Y >=0.263, with a la 
spread <x(Y)=0.010. The whole HB is characterized by 
< Y >=0.265, with a la spread ct(Y)=0.010. This means 
that the full range of Y values is already present in the 
RHB sample. In the reasonable assumption that the spread 
of initial Y scales with the extension of the Na-0 anti- 
correlation, we should expect to find among RHB stars the 
full range of Na-O abundances observed along the RGB. 
A range of He abundances is necessary along the RHB to 
match the observed thickness in m,F33&w, and the large 
range of AM (at fixed Y) is required to match the ex- 
tended and almost horizontal structure of the RHB in the 

(m_F336lV, WF336 — WF555W ) CMD. 

Figure [13] displays the mass and Y distribution as a 
function of the (m,F336w — fiFsssw) colour in one represen- 
tative synthetic realization of the cluster HB, compared to 
the M13 counterpart. The range of Y and the shape of the 
trend with colour are very different. The maximum mass of 
HB stars in M79 reaches higher values compared to M13, 
and the minimum HB mass is also higher than in M13. 
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Table 1. Filters adopted, number of observed HB and RRLyrae 
stars and values of Y max and < Y > for the three clusters of our 
analysis. 



Lyme 



< Y > 



M3 F336W 158 

F555W 
M13 F160BW 228 

F336W 

F555W 
M79 F160BW 130 

F336W 

F555W 



59 0.266 
7 0.300 

5 0.281 



0.256 
0.270 

0.265 



6 SUMMARY AND CONCLUSIONS 

There is a general consensus that the first parameter shaping 
the globular cluster HB morphology is metallicity and the 
second one has been suggested by several authors (Lee et al. 
2002; G10; D10) to be the age. While in general this picture 
may be surely adequate, however it is able to explain only 
few young clusters in the Galaxy, while it fails to reproduce 
a number of "classical" cases where an additional ingredient 
(sometime called third parameter) is required. The "third 
parameter" has been proposed to be the stellar density or 
luminosity density (Fusi Pecci et al. 1993; D10) or the initial 
helium abundance (G10). 

In this framework we considered three GGCs with similar 
age and metallicity, M3, M13 and M79, with the aim of 
studying their HBs. According to G10, M3 is younger by 
~ 2Gyr than M13, however we find, in agreement with a 
number of authors (Salaris & Weiss 2002; Marin-Franch et 
al. 2009; D10), that these clusters differ in age by less than 
1 Gyr. 

This triplet is well known to have extremely different 
HBs, with M13 displaying by far the bluest morphology, 
and showing also evidence of at least 2 gaps at different 
effective temperatures, while M3 has been proposed over 
the years as a prototype of a "normal" HB. 
As done in Paper I, we compared HST UV and optical 
CMDs with theoretical HB models, with the aim of un- 
derstanding the impact (if any) of cluster-to-cluster Y 
variations on HB morphologies. At odds with the case of 
NGC2808, the three GCs analyzed in this work, do not 
show (so far) evidences of quantized MS, so our analysis 
has been performed without any a priori knowledge of the 
Y distribution. As explained in Section 4, this forced us to 
enlarge the parameters space: the quantities needed are the 
minimum value of Y (Y m i n ) and its full range (AY), the 
mean value of the mass lost along the RGB AM and the 
dispersion around the mean value. As a consequence, it is 
hard in these cases to determine uniquely a full and detailed 
representation of the Y variations along the observed HBs. 
A more stringent and robust derivation involves instead the 
largest value of Y (Y max ) that is clearly constrained by the 
distribution of stars in the (mFisoBW', mFWOBW — mF555w) 
CMDs. For the case of M3, given the shorter extension 
of its HB, the magnitudes mp336w give already a strong 
constraint on Y max . 
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Figure 13. Mass and Y distribution as a function of the 
( m F336W — m F555w) colour along the HB of M79 (grey crosses) 
and M13 (black dots). 



With these caveat in mind, we performed a detailed 
analysis by comparing observed CMDs with a large set 
of synthetic CMDs and ZAHBs from the BaSTI database 
(Pietrinferni et al. 2006). The effects of levitation for 
stars with T e g > 12000 K and anti-correlations along the 
HBs (Marino et al. 2011; Gratton et al. 2011) have been 
considered. 

We find differences AY max ~ 0.02 — 0.04 between these 
GGCs. In particular M13 displays the largest value 
(Y max ~ 0.30), M3 (Y max ~ 0.27) the smallest one, and 
M79 is an intermediate case with Y max ~ 0.28. These values 
are summarized in Table 1. They seem to qualitatively 
correlate with the differences in the temperature (colour) 
extensions of the cluster HBs. It is also interesting to 
note that our estimates of Y max for these three clusters 
nicely correlate with the observed range of light-element 
variations. In particular, M13 shows the most extreme 
Na-0 anti-correlation (Sneden et al. 2004), with stars 
reaching [0/Fe]= -1.1 and [Na/Fe]= 0.7, while M3 the 
least extended one, [0/Fe]= -0.2 and [Na/Fe]= 0.5. This is 
in line with the strict correspondence between HB colours 
and Na-O abundances observed in NGC2808 by Gratton et 
al. (2011). 

The comparison between M13 and M3 is particularly 
interesting, because of several previous analyses of their 
HB morphology and Y distribution. For example Catelan 
(2009) performed a detailed analysis comparing the mean 
masses of the HBs of these clusters with age differences 
proposed by different authors. By using several mass-loss 
recipes, he found that there is no way to reproduce the 
different HBs of M3 and M13 only in terms of age, but at 
least one additional ingredient should be invoked to account 
for the blue HB extension observed in M13. Our best-fit 
model for M13 requires to split the simulation in three 
steps, given that no single continuous distribution of Y and 
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mass loss allows to match the observed CMDs. Thus we 
divided the HB in three groups: 1) the RHB with a typical 
Y min = 0.246 and a mean mass lost = AM = 0.21M Q , 

2) stars with —3 < {mF\§QBw — TOfsssvk) < —1.5 best 
reproduced with < Y > = 0.285 and AM = 0.235M© and 

3) stars with (rriFieoBw — rnFB&sw) < —3 which show a 
< Y > = 0.30 and AM = O.266M . A fit with a population 
of stars with uniform Y — 0.265 cannot be ruled out on 
the basis of only the HB analysis, while a single population 
of stars with Y = 0.28, as proposed by Caloi & D'Antona 
(2005) and D'Antona & Caloi (2008), is incompatible with 
the distribution of stars in our UV CMD. 

For M3 a single synthetic population with Y m i n = 0.246 
and distributed according to a uniform probability distri- 
bution with AY = 0.02 is required. A total mass loss of 
AM = O.122M0 and a linear increase as a function of Y, 
as constrained by the fit of star counts as a function of 
magnitude and colours, is needed. In the Appendix B we 
will discuss the Y distribution derived from our simulations 
in the context of the constraints posed by the analyses by 
Catelan et al. (2009) and D'Antona & Caloi (2008); it turns 
out that our estimate of Y max = 0.266 is robust. 
As highlighted by G10, while M13 seems to behave as other 
relatively massive clusters, M3 appears to be peculiar and 
a more extended HB would have been expected in this 
case. The age difference proposed by G10 between these 
two clusters would justify at least in part the different HB 
morphology, in the framework in which age is the second 
parameter. Our results would instead lead to think that M3 
and M13 experienced a different amount of enrichment of 
light elements. This would be compatible with the scenario 
proposed by G10 (see also Carretta et al. 2009) that invokes 
a delayed cooling flow in the case of M3. In particular the 
HB simulations and derived Y distributions would suggest 
that M13 is qualitatively similar to NGC2808, and that it 
likely experienced a similar star formation, while M3 (and 
M79) probably had a less complex formation history. 
Analyses based on a suitable combination of UV to optical 
photometry - and when available observed RR Lyrae period 
distributions - and synthetic HB simulations can provide 
not only insights on the HB second parameter problem, but 
in general they can potentially give some clues on the first 
stages of GCs formation and chemical evolution. 
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Appendix A 

The effect of anti-correlations on bolometric 
corrections 

We checked the impact of CNONa anti-correlations on 
the a-enanched bolometric corrections by performing a num- 
ber of tests. In particular, we have chosen three represen- 
tative pairs of (T e fi, Iog(g)) values along our theoretical 
Q-enhanced ZAHB with [Fe/H]=-1.62, Y=0.246, namely 
(5600 K, 2.40), (7000 K, 2.80), (10000 K, 3.50). For each 
of these pairs we have calculated a synthetic spectrum with 
the SYNTHE code by R.L. Kurucz (Sbordone et al. 2005) 
in the wavelength range between 1200 and 7500 A, includ- 
ing all the atomic and molecular lines available in the lat- 
est version of the Kurucz/Castelli linelist. The spectra were 
computed employing the LTE, plane-parallel model atmo- 
spheres calculated with the ATLAS12 code (Castelli 2005). 
ATLAS 12 adopts the "opacity sampling" method in the 
line opacity calculations, and allows to generate models for 
arbitrary chemical abundance mixtures. All the model at- 
mospheres and synthetic integrated spectra were computed 
with [Fe/H]=-1.62 and Y=0.246 and both the standard 
[a/Fe]=0.4 metal distribution used in the stellar evolution 
models, and the CNONa2 mixture by Sbordone et al. (2011). 
This latter mixture displays - compared to our reference a- 
enhanced mixture - enhancements of N and Na by 1.44 dex 
and 0.8 dex by mass, respectively, together with depletions 
of C and O by 0.6 dex and 0.8 dex, respectively, and un- 
changed CNO sum. This mixture is representative of ex- 
treme values of the CN and NaO anti-correlations observed 
in globular cluster stars (see Sbordone et al. 2011 for de- 
tails). It is also very important to notice that low- mass stel- 
lar evolution models calculated for this mixture turn out 
to be identical to models calculated with the standard ot- 
enhanced metal distribution. 

From the synthetic spectra we have calculated bolo- 
metric corrections to the F160BW, F336W and F555W fil- 
ters used in our analysis, for both metal mixtures. Only 
at T cff =5600 K there is a difference for the F336W fil- 
ter (we did not consider the F160BW at this low tempera- 
ture, because our photometry does not register cold HB stars 
in this passband) that amounts however to only 0.03 mag 
(CNONa2 mixture stars appearing fainter). Given the small 
value of this difference and the fact that such extreme val- 
ues of the anti-correlations appear not to be reached at such 
low temperatures along the HB, the outcome of our analysis 
of NGC2808 and the results of this paper, based on stan- 
dard Q-enhanced bolometric corrections, are unaffected by 
the presence of the observed CNONa anti-correlations. 

We have explicitly determined also the effect of He 
variations (at fixed metal content) on the bolometric 
corrections. We considered the point with (7000 K, 2.80) - 
i.e. around the blue edge of the RR Lyrae instability strip - 
and calculated spectra for our standard a-enhanced metal 
distribution, with Y=0.35. Compared to the case with 
standard Y=0.246, both F336W and F555W bolometric 
corrections are changed by less than 0.01 mag. The F160BW 
filter is affected at this level by less than 0.02 mag. 
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Figure 14. Comparison of M3 HB with a ZAHB sequence with 
Y=0.256, for log(T eff ) > 3.89 (see text for details). 



Appendix B 

Comparison with previous studies on the HB 
morphology of M3. 

We address here two issues that arise from a comparison 
with recent results from the literature about the HB mor- 
phology of M3. First of all, the range of Y derived from our 
simulations is about twice the value determined by Cate- 
lan et al. (2009) from a comparison of theoretical ZAHB 
magnitudes and gravities with both Stromgren photometry, 
and spectroscopically measured surface gravities and T e s. 
The size of the sample with spectroscopic gravities and T e g 
estimates is relatively small, and within the corresponding 
error bars models with Y ~ 0.27 cannot be easily ruled 
out. More compelling seems to be the constraint from the 
Stromgren CMDs. A fit of theoretical ZAHB models to the 
magnitude (y, b and v) of the lower envelope of the observed 
HB at both the red side of the instability strip and at the 
'knee' of the stellar distribution beyond the blue edge of the 
strip, before the sparsely populated and not very extended 
blue vertical tail - where ZAHBs of different Y tend to be- 
come degenerate - discloses a range of Y at most equal to 
AY^O.01, with Y increasing at the blue side of the strip. 
This disagreement with our results is however only appar- 
ent, for it can be very easily explained by recalling that to 
any given Y value corresponds it is associated a dispersion 
in colour along the HB (see Section 4.1). Figure [TJ] illus- 
trates clearly this point by comparing a theoretical ZAHB 
for Y=0.256 and T e fi larger than the blue edge of the insta- 
bility strip - shifted by our derived distance modulus and 
assumed reddening - to the observed data. Notice how for 
(0.2 < (niF336w — mF555w) < 0.4), that corresponds to the 
region of the 'knee' in the Stromgren CMDs, this ZAHB fits 
nicely the observed lower envelope of the CMD, even though 
the average Y in that colour range inferred by the best-fit 



synthetic model, is 0.260. Given the existence of a range of 
HB masses at fixed Y, this colour interval is also populated 
by stars with Y as low as ~0.255, that are on average the 
faintest objects and hence determine the observed lower en- 
velope of the magnitude distribution. 

The second issue is related to the cluster rich 
RRLyrae population. The period distribution of fundamen- 
tal RRLyrae is strongly peaked, and in the recent literature 
attempts have been made to reproduce both the peaked pe- 
riod distribution and the HB colour extension. Castellani, 
Castellani, & Cassisi (2005) have shown that when keeping 
Y constant, a suitable bimodal mass distribution with two 
different mass dispersions is able to reproduce the observed 
features. On the other hand, D'Antona & Caloi (2008) as- 
sumed a range of Y and a very small mass dispersion for 
the mass loss along the RGB, and reproduced the observed 
HB colour extension and RRLyrae periods by calibrating 
the Y distribution (instead of the HB mass distribution). It 
is interesting to notice that D'Antona & Caloi (2008) find 
from their simulations that 90% of HB stars are expected 
to have Y^JO.265, in agreement with our results, and only 
10% of objects in their simulations reach Y values above our 
estimated maximum Y. 

Here we have considered the synthetic objects in the 
best fit simulations that fall within the boundary of the in- 
stability strip (using the same T e fi boundaries employed by 
D'Antona & Caloi 2008) and were excluded from the com- 
parison with the observed CMD. We have calculated their 
pulsational periods using the relationships by Di Criscienzo 
et al. (2004) and compared with the observed distribution 
(Corwin & Carney 2001). Our predicted period distribution 
does not reproduce the observed ones, which is probably not 
surprising given our 'simple' assumptions about the distri- 
bution of Y among stars fed onto the HB. 

As a test, we have calculated additional synthetic HB 
models for M3, employing this time the adhoc Y distribu- 
tiorjf] proposed by D'Antona & Caloi (2008 - see their Fig- 
ure 4) to reproduce both the observed period distribution in 
the instability strip and the (B—V) colour distribution along 
the HB, while keeping the same mass loss for all HB stars. As 
in D'Antona & Caloi (2008) we employed a Gaussian mass 
distribution; our adopted mean mass loss is AM=O.16OM0 
with a small (as discussed in D'Antona & Caloi 2008) Gaus- 
sian spread cr( AM)=0.002Mq , and employed the same red- 
dening and distance modulus as in the previous simulations. 
These choices allow to reproduce well the narrow period dis- 
tribution (both mean value of the logarithm of the periods 
and the la dispersion) whereas a comparison of the the stel- 
lar distribution in the CMD at the blue side of the strip is 
less satisfactory. The left and right panels of Figure [TJ] show 
the observed and simulated (iriF336w, rn,F336w — rriFsssw) 
CMD. Notice how the observed shape of the CMD is rea- 
sonably reproduced as long as Y is below ~0.270. For higher 
values of Y, synthetic stars tend to reach too bright magni- 
tudes (while in the BV plane they tend to be located along 
the more vertical blue part of the HB, at/beyond the 'knee'). 
This reinforces our estimate of the Y range as determined by 
our previous set of simulations - that is the main purpose of 
our analysis - and the need to increase the mass loss spread 



Our minimum Y value is 0.246 instead of 0.240 
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Figure 15. Observed (left panel) vs one realization of a syn- 
thetic (right panel) (m F33ew , m F33ew — m F555w ) CMD of M3 
using the Y distribution proposed by D'Antona & Caloi (2008). 
Open symbols denote synthetic objects with Y>0.270. The ZAHB 
sequences for Y=0.246 and Y=0.300 are also displayed. The rect- 
angular box marks the region considered in the comparison (see 
text for details). 

at the blue side of the instability strip - instead of Y - in 
order to match the CMD location of the bluer objects. 
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